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>tis shown that for collisionless transverse magnetosonic shocks, there exists an intermediate 
range of Mach numbers between the maximum Mach number for transverse resistivity and the 
onset of significant ion reflection. In this intermediate range, the structure is dominated by a 
breakdown of the quasi-neutrality approximation and longitudinal resistivity is a simple fluid 
model for the dissipation mechanism which forms the shock. 
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LONGITUDINAL RESISTIVITY AND INTERMEDIATE MACH NUMBER 
COLLISIONLESS TRANSVERSE MAGNETOSONIC SHOCKS 

I. Introduction 

It is now reasonably well established that transverse collisionless 
shocks can exist in at least two Mach number regimes. For m (the 
magnetosonic Mach number) less than about 2.5, transverse resistivity is 
the dissipation mechanism which forms the shock.Since the plasma is 
nearly collisionless, this resistivity is anomalous in that it is assumed 
to arise from an electron-ion streaming instability in the shock front. 
Measurements of the turbulent spectrum in the shock front show that 
fluctuations with wavenumber k ~ kp, the Debye wavenumber, are indeed 
excited.It is generally assumed that the resistive dissipation heats 
only the electrons. However a small amount of ion heating cannot be ruled 
out. It is certainly not easy to see how an instability can heat only one 
species. An ion acoustic instability, at least, does give rise to a small 

amount of ion heating. ^ In fact one theory of the resistive shock 
o 

structure 0 has this small ion heating playing an important role in 
determining shock structure. 

For larger m, it is now well established both from theory^’and 
particle simulations,12-19 that ion reflection from the shock front can act 
as a dissipation mechanism in one of two ways. On the one hand, these 
reflected ions may stream ahead of the shock if they are unmagnetized or 
sufficiently weakly magnetized.9,10,12-15 If thig be the case> the 
standard kankine-Hugoniot conservation relations do not apply because of 
the additional species upstream. On the other hand, the reflected ions may 
turn around in the magnetic field and drift downstream, their large Larmor 
orbit velocity effectively becoming part of the downstream 
pressure. 11 ' 1 ** -1 ^ Since there is no additional upstream species in this 








case, Rankine-Hugoniot conservation relations do apply. In fact the 
simulations of Refs. 18-19 were successful because the Rankine-Hugoniot 
relations were exploited in initializing the system. Of course for the 
latter dissipation mechanism to be significant, the time of the experiment 
or simulation must be longer than about ir/u) cl , where u) ci is the ion 
cyclotron frequency. 

The purpose of this paper is to show that between these two ranges of 
Mach number, there exists an intermediate range. In this intermediate 
range, the shock structure and dissipation are dominated by the breakdown 
of the quasi-neutral approximation. Above a critical Mach number, the 
quasi-neutral fluid equations predict infinite slope, implying, of course, 
the breakdown of quasi-neutrality. However, the full set of fluid 
equations, including charge separation, are nonsingular. Thus, above this 
Mach number, the flow sets up an electrostatic ion acoustic wave train. 
The Landau damping of this wave is the dissipation mechanism which forms 
the shock. This Landau damping may heat either the electrons or ions. In 
purely electrostatic shocks with very cold ions, the electron Landau 
damping seems to dominate. 20 Since the electron Landau damping rate is 
small, the shock structure is oscilliatory. On the other hand, in 
magnetosonic shocks the ions are warmer, so ion Landau damping dominates 
and the shock structure is monotonic. It is interesting that this 
dissipation mechanism does not rely on exciting some plasma instability in 
the shock front, and in that sense it is not anomalous. Indeed, the 
electrostatic subshock has a Debye length scale size, so it is almost 
impossible that an instability with still smaller wavelength could be 


excited. 








While this dissipation mechanism is not anomalous it is very 
complicated* The simplest model of it shows tnat it is nonlocal, and in 
that sense non fluid like. However, there is a new simple fluid 
dissipation model. When the electrostatic oscillation is set up, the 
longitudinal electron and ion velocities are no longer equal. A 
longitudinal friction force between the two is an additional possible 
dissipation mechanism which apparently has not been examined in the 
previous literature on high Mach shock waves. An important purpose of this 
paper is to investigate the effect of such longitudinal resistivity on 
shock structures. Our conclusion is that it is a very good fluid model for 
collisionless magnetosonic shocks with Mach number in the intermediate 
regime. It has a precise onset, namely when the shock reaches the critical 
Mach number. It shows that shocks above this Mach number have an 
electrostatic subshock embedded within in and its structure is governed by 
longitudinal resistivity. When the Mach number increases to a second 
critical value, the full fluid equations become singular, and this defines 
the upper limit of the intermediate range. The interpretation of this 
second singularity is simply made in terms of ion reflection. Above this 
point, the shock, if it exists, is governed by ion reflection as described 
in, for instance. Refs. 9-19. At no point is it either necessary or 
physically correct to postulate a viscosity. 

While the theory we deal with here concerns only steady state shocks, 
the general conclusions should also apply to non-steady flow of 
interpenetrating plasmas. Specifically, electrostatic interactions should 
be important in a recent series of laboratory experiments on the flow of a 
laser produced plasma through a background magnetized plasma. 









In Section II we review experimental data 


high Mach shock 


structures and show that It Indicates the existence of this intermediate 

regime* To do so, we note that in this Mach number range many of the 

experiments are over in a time of order of or less than tt/u> , , so 

ci 

1 3 22 23 

convection of reflected ions downstream * »» cannot be the dissipation 
mechanism of importance* Also, downstream ion temperatures are measured in 
several experiments and they agree with Rankine-Hugoniot. This means that 
pure reflection of unmagnetized ions,^*^ cannot be the dissipation 
mechanism either. 

In Section III we discuss shock structure for unmagnetized ions where 
the quasi-neutrality assumption is valid and also where it is not. V7e find 
that for sufficiently weakly magnetized plasmas, characteristic of the 
experiments in Refs. 1-6, 22, 23, the electrostatic and magnetic structures 
essentially decouple, and that the electrostatic subshock must ride on the 
crest of the magnetic shock structure. This also agrees with experimental 
results. 

In Section IV, and also In Appendix B, we discuss the actual 
dissipation mechanism (Landau damping) and fluid models for it. As we 
show, even the simplest approximation to it is quite complex, and does not 
easily fit into a simple fluid model. However the actual physical model 
provides Important insight on a best guess for longitudinal resistivity and 
on the partition of resistive heating between electrons and ions. Another 
interesting phenomena is that Landau damping can only be a valid long time 
description of the dissipation mechanism if the shock structure fluctuates 
in time. Thus the shock cannot be time independent but, of course, its 
average profile can be. 
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Section V presents the results of shock structure calculations in the 
intermediate Mach number regime. Also, where possible, comparisons are 
made with experimental results. Finally, Section VI draws conclusions. 








II. Review of the Results of Laboratory Experiments 


on Transverse Collisionless Shocks 


In this Section we review briefly previous experiments of transverse 
shocks in laboratory plasmas. 1 “6,22,23 our conclusion, based on 
examination of this data, is that for m < 2.5, the shocks are resistive, as 
described by authors of these references. However, there is a 
range 2.5 < m < 5 where steady shocks form which are above the critical 
Mach number for resistive shocks. However, the process of ion reflection 
followed by downstream convection by the E_ x B_ drift cannot be operative 
either, if only because the duration of the experiment is less than the ion 
Larmor time. Furthermore, the dissipation mechanism cannot be simple ion 
reflection of magnetized ions either. For ion reflection in an 
unmagnetized plasma, the second ion species upstream drastically alters the 
Rankine-Hugoniot conditions downstream. 10 However, in two of the 
experiments cited,d owns tream ion temperatures were measured and found 
to agree well with Rankine-Hugoniot conditions. Thus our conclusion is 
that there must be some sort of new dissipation mechanism associated with 
this intermediate range of Mach numbers. 

In Ref. 1, Paul et al . launch shocks with Mach numbers 2.5, 3.7 and 
6.3. For the first, resistivity forms the shock and gives the electron 
heating required to satisfy the Rankine-Hugoniot relations. In the case 
of m • 3.7, the downstream electron temperature is less than that predicted 
by Rankine-Hugoniot, implying that the ions are heated. In addition to the 
normal magnetic structure, there is a very low amplitude magnetic pedestal 
in front of the shock. This pedestal is generally interpreted to mean a 
small number of ions are reflected from the shock front and convect the 
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magnetic field upstream. However in Paul's experiment, the convection of 
these ions downstream by the E_ x B_ drift is almost certainly not the 
mechanism which forms the shock at m * 3.7. The shock velocity is about 
2.5 x 10 7 cm/sec and the shock propagates ahead of the piston a distance of 

10-15 cm, implying a time of about 400 nsec. However, ir/u)^ is also about 

400 n sec. Thus any reflected ions cannot make it back behind the shock 
until the shock hits the cylindrical axis. However, shocks are observed 
for much shorter times. For the measurement at M - 6.3, the upstream 

pedestal broadens and increases in amplitude so that it encompasses the 
entire shock. At this Mach number, reflected ions are almost certainly 
playing a role. 

In another series of experiments, Keilhacker et al.^ measured the 

properties of magnetosonic shock waves with Mach numbers between 1.5 and 

5. Here ion temperatures behind the shock are measured by collective 

scattering of laser light at wavelengths less than the Debye length. The 
measured ion temperatures agree with Rankine-llugoniot conditions. In this 
experiment, the total transit time for the shock is about ly sec. Thus the 
dissipation mechanism cannot be the convection of reflected ions downstream 
because there is not enough time; neither can it be the pure reflection of 
unmagnetized ions, because the downstream ion heating corresponds to that 
predicted by Rankine-Hugoniot. 

Segre and Martone^ have also investigated transverse shocks in the 
Mach number range 3.1 < m < 4.25. In their experiment, there is no 

magnetic field pedestal leading the shock, so reflected ions almost 

certainly play no role. Even if there were reflected ions, the ion period 
ir/u c ^ ~ 450 n sec, the propagation time is only about 250 n sec, so no 
reflected ions could make it back downstream. In Ref. 22 the downstream 


7 














electron density and temperature are measured by laser scattering. However 
magnetic measurements are used to measure n(T g + T^) downstream, so that 
downstream ion temperatures can be directly inferred. The measured ion 
temperatures do in fact agree with that predicted by Rankine-Hugoniot. 

Another series of shock experiments are those of Eselevich et al.^ 
The shocks are produced in a tube of 8 cm radius in H, He or A. Since the 
gas can be varied, and the density varies between 10^ < n < 10^ cm ^ and 
100 G < B < 2000 G, there was a very large parameter space studied. 
Eselevich et al. found that for Mach numbers between 3 and 5, steady she s 
form. In this range of Mach numbers, a pedestal in magnetic field might 
might not lead the shock. For instance, Fig. 2a of Eselevich et al. si 
a Mach 4 shock in helium with no pedestal, whereas Fig. 4 shows a Mach 4 
shock in hydrogen with a small (i.e., less distinct than in Paul's^ - 
experiment) magnetic foot leading the shock. Thus, reflected ions are not 
a universal, and probably not a very important, aspect of steady transverse 
shocks in this intermediate range of Mach numbers. However, Eselevich et 
al., like Paul et al., did observe that as the Mach number increases 
(beyond about 5.5 generally in Ref. 23), the magnetic foot does expand out 
and dominate the entire structure. 

The most distinct feature of Eselevich's observations are that in the 
intermediate range of Mach numbers, an isomagnetic potential jump forms on 
the crest of the shock wave. The thickness of the potential jump is 
typically a Debye length. Thus, a key feature in this intermediate range 
is the breakdown of the assumption of quasi-neutrality. The potential jump 
across this isomagnetic transition region initially increases with Mach 
number, but then above some higher Mach number begins to decrease until 
above a second critical Mach number, typically about 5 or 6, the entire 






isomagnetic jump disappears and the shock structure is taken over by the 
magnetic pedestal. Hence, we feel that a key element in a fluid 
description of shocks in this intermediate range is the breakdown of quasi- 
neutrality. Once quasi-neutrality is violated, an additional fluid 
dissipation mechanism, longitudinal resistivity becomes possible. In the 
next two sections we examine the consequence of longitudinal resistivity. 
Although the actual dissipation mechanism is considerably more complicated 
and non fluid like (see Sec. IV), a longitudinal resistivity model can 
explain many features of the experimental results. Specifically it hinges 
on the breakdown of quasi-neutrality; it provides for the proper electron 
and/or ion heating; it predicts an isomagnetic potential jump on the crest 
of the shock wave; and it gives the shock structure only in an intermediate 
range of Mach numbers between the critical Mach number for transverse 
resistivity and the onset of strong ion reflection. Thus, longitudinal 
resistivity is probably as good an approximation to the dissipation 
mechanism that one can obtain within the limits of a conventional fluid 


description. 










In this Section, we discuss shocks formed by longitudinal resistivity. 
Before doing so, however, in Section III.A we set the stage by reviewing 


the structure of shocks with transverse resistivity. Section III.B 
discusses the effect of longitudinal resistivity, and Section III.C sets 
out the procedure for calculating the shock structure in weakly magnetized 
plasmas. 
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A. Review of Shocks Formed by Transverse Resistivity 

To start we formulate the equation for quasi-neutral transverse 
resistive shocks. We neglect both transverse and longitudinal electron 
inertia, and the motion of the ions in the transverse direction. Assuming 
that the shock propagates in the x direction, the magnetic field is in the 
z direction, and upstream parameters have a zero subscript, then the 
electron and ion density equation gives the result 


( 1 ) 


where a velocity or electric field without a direction subscript means a 
velocity or electric field in the x direction. Then the electron momentum 
equation in the x and y direction give respectively equations for E and 
V ey . They are 


E 


1_ d_ 
ne dx 


n 


T 

e 



( 2 ) 


and 


eV n 

(b -_°I 

mvc K o n 


(3) 


where m is the electron mass, v the effective collision frequency, and the 
resistive force density on the electrons is - nmv(V ey - ) = - nmv V gy 

since we assume V » 0 (unmagnetized ions). Also we have used the fact 


To write out the electron temperature equation we first assume that 
the electrons have y - 5/3. Then we must also specify what fraction of the 
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resistive heating goes into the electrons and what fraction goes into the 
ions. In terms of the ion y, the total energy equation for electrons plus 


fi nMV 3 + 1 nVT + -X. nVT. ) - - neV E , 
dx '•Z 2 e y-1 L J ey y’ 


whereas the total energy equations for the electrons and ions separately 


(1 nMV 3 + ^ nVTj - neVE - A, 


where M is the ion mass. The term A denotes an energy exchange between 
electrons and ions, and its value depends on what fraction of the resistive 
heating goes to electrons and what fraction goes to ions. To determine 
this, one must write out the temperature equations for electrons and 
ions. To do so, take the dot product of the electron momentum equation 
with Vg and subtract the result from the electron energy equation to get 


- nVT - V 4 - nT - A + nmvV 2 
e dx e ey 


Then multiply the x component of the ion momentum equation 


^ (nMV 2 + nT i ) - neE 


(7) 









by V, subtract the results from the ion energy equation to obtain 


T lx nVT i ” v nT i " ” A ' 


If a is the ratio of ion resistive heating to electron resistive heating. 


A - - nmvV 2 . 

1 + a ey 


The quantity a depends on the physical process causing the resistive 
heating. Classical resistivity gives a - m/M. Anomalous resistivity gives 
an a which depends upon the particular mechanism. Generally <* « 1; for 
Instance if the resistivity is caused by an ion acoustic instability, 
a - 0.1-0.2.^»® Finally, these equations are closed by Ampere's law 


In order to see where the critical Mach number is, a standard 
procedure is to linearize Eqs. (1-10) about the upstream or downstream 
state. Assuming that ~ ■ 1c, and an unsubscripted variable denotes either 
the upstream or downstream state, there are only two roots for k, 
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k - 0 


(11a) 



k > 0 so the linearized solution diverges away from the state. Equation 
(12) corresponds to the upstream state of the shock wave for which the flow 
speed must exceed the magnetosonic velocity. The second regime corresponds 
to 


where k < 0. Since the flow speed is less than the magnetosonic speed, it 
corresponds to the downstream state, and since k < 0, perturbations 
converge into the downstream state. This then is a standard (transverse) 
resistive shock having a subcritical Mach number. The final regime 
corresponds to 
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( 14 ) 


having k > 0. Since perturbations now diverge away from the downstream 
state, there can be no shock in the regime defined by Eq. (14) if 
transverse resistivity is the only dissipation mechanism. This then 
defines the critical point above which no shock can form. 

To continue, let us write out the equations for a quasi-neutral shock 
as a system of nonlinear first order differential equations, as is 
appropriate for standard numerical integration. 

The equations for B and V ey are already expressed in the proper format 
in Eqs. (10) and (3). Equations (6a) and (6b) can easily be reduced to 


but now dn/dx appears on the right hand side. To solve for dn/dx, use the 

ion momentum equation, insert the electric field from Eq. (2), use Eqs. 
dT E dT 

(15) and (16) for ^— and t— , and solve for -r— . The result is 


2 2 

V - nmvV 

: -££b+|—^ 

c 3 n V 

_ o o 

m 2 - yt - 1 t 
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This can be used on the right hand side of Eqs. (15) and (16). Notice 
. dT. dT 

that dx -an< ^ E become infinite at V - (yT^ + -y T e )/M. Thus, 

2 

quasi-neutral shocks formed by transverse resistivity only exist if V < 
(yT^ + y T^J/M everywhere. If the shock is monotonic and y * -y (or more 
generally, if y g * y^), this critical Mach number can be determined by the 
Rankine-Hugoniot conditions only. However, if the shock is not monotonic, 
and/or y^ * y^, the critical Mach number depends on the details of the 
oscillatory structure and how the resistive heating is partitioned between 
electrons and ions. 

A key question then, for the description of shocks in the intermediate 
regime, is the physical significance of the singularity. One might think 
that the singularity indicates the onset of ion reflection, but this is not 
the case. To see this, imagine that T^ « T 0 . Then ion reflection occurs 

where V ~ (T, /M)*^ - 0. However, the singularity occurs at a much higher 
5T e 1/2 

velocity, V • (y^—) . The significance of the singularity lies not in 
ion reflection, but rather in the breakdown of the quasi-neutral 
assumption. This is clear from the fact that both E and , as expressed 
in Eq. (3), become infinite at the singularity, implying a breakdown of 
quasi-neutrality. Thus, in the intermediate regime, charge separation 


begins to play an Important role in the shock structure. 
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B. Shocks with Longitudinal Resistivity 

In this subsection and the next, we discuss how shocks are formed by 
longitudinal resistivity. If quasi-neutrality is no longer assumed, the 
steady state ion fluid equations are 


37 (“i^i + n i T i) " n i eE “ n i m M V i - V* (19 > 

and 

h (7 “i m i + -iVi) ■ n iV E + T+a V<r 

< 2I » 

where and are defined analogously to a and v. Notice that we allow 
different values for longitudinal and transverse dissipation and energy 
partition. This is reasonable because, as we will see, these two types of 
dissipation arise from completely different physical processes. 


The fluid equations for the electrons are 







■ — n.mv fl fv, - v ) 2 . 

1 + a„ i i e J 


V --- (B - — B). 

ey mvc o n 1 


These are supplemented by Maxwell's equations 


(23) 


(24) 


4im 

i® -_® v 

dx c ey 


(25) 


and 


|| - 4 ir e(n 1 - n @ ). (26) 

As in the previous section, we begin by linearizing these equations 

about the upstream or downstream state^ As before, there^is no resistive 
2 ”e "i 

heating in linear theory, so - -j T g — and - (y-lJT^ — . Of course, 
^ * - V n @ /n and ^ - V n^/n. Then linearizing, a straightforward 


calculation gives the result 









k 2 ( YTi -MV 2 )(f k 2 T e 


a) V 

--* 5 - 2 . 


^ v # kV ? c mV u) 2 . , 

- ‘"» 0 « (1 + - V 2 -) 11 It T e + ^ T 1 - -V s —T— + n T i - av ] ■ 

w pe 


k-^v- 


(27) 


where w c is the electron cyclotron frequency. Factoring out the two roots 
at k • 0, what remains in Eq. (27) is a cubic equation for k. As we will 

2 2 

see, one of the roots is of order ui V/vc , while the other two roots are 

pe 

2 1/2 2 2 
of order k D = (4^/1^) ' . Assuming that k D » oi^V/vc , the equation 

separates. The low k root is still given by Eq. (lib). We denote this the 

Magnetic Mode. To get a simple form for the high k roots, we further 
2 2 

assume that muj MkpV V q « 1, and <i) c /v « 1. Since the length of a 
resistive shock 1-5 is typically 10c/m pe , an approximate order of magnitude 
for the anomalous collision frequency is v - 10(D pe V/c. The first 
inequality is then easily satisfied for the parameters in the experiments 
of Refs. (1-6, 22, 23). As for the second inequality, an estimate is 
required for v^. As we will see in the next section, is very large, 

larger even than u) pe » Thus the second inequality is also satisfied. In 


this case, the solution for k becomes 









tCt*.**v 


¥■ (** - 


™ [^-MV 


■ l " TT K-T-[4 T e + * T i - M ^l) 2 + 


T^tf T « + K T 1 - MV ) (t t i - ™J)| l 


where u £ is the electron thermal velocity u g » (T e /m)^^. 

Notice that the wave number given by Eq. (29) is of order k D . Also, 
the magnetic field completely decouples from this large k mode, so it is 
entirely electrostatic in nature. Depending on the values of velocity and 
temperature, the linearized modes structure themselves as shown in Table I. 

Now consider how to integrate Eqs. (18)-(26) to solve for the shock 
structures. Let us say we start upstream and integrate toward the 
downstream state. Referring to Table I the upstream plasma is in Class A, 
since the flow speed must exceed the magnetosonic speed. Hence there are 
two linearized perturbations which increase away from the upstream state. 
The relative amplitude between the two is arbitrary at this point. Now let 
us assume the downstream plasma is in Class B or C. In each case there are 
two linearized eigenfunctions which converge into the downstream state, and 
one which diverges. Thus the relative amplitude of the two linearized 
perturbations upstream must be chosen so as to eliminate the divergent 
perturbation downstream. This is then a nonlinear eigenvalue problem, and 
this concept has been applied to both conventional shocks^ and reflected 
10 


ion shocks. 








Table I 


Velocity Regimes and linearized modes in each regime. 


Velocity Range 


+ T + yT 1 < ! 


'"'A + T T e + * T 1 


C yT i < MV 2 < T e + 3T i 


D MV < yT 


One electrostatic mode and the 


magnetic mode diverge from the 


increases in the 


direction of flow. 


electrostatic mode converges. 


Magnetic mode converges, one 
electrostatic mode converges, one 


Magnetic mode diverges, both 
electrostatic modes converge. 


Magnetic mode and one electric mode 
diverge, the other electrostatic 


mode converges. 
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If the downstream state is Class B, the scheme is trivial. The 
divergent downstream perturbation is electrostatic so the solution is just 
to eliminate electrostatics and calculate shock structure based on quasi¬ 
neutrality. If the downstream state is Class C, the problem is more 
complex. Here the divergent state is magnetic and both convergent states 
are electrostatic, so there is necessarily a coupling between the magnetic 
and electric structure somewhere in the shock. In the next subsection we 
will show how to calculate this coupled structure in a simple way by 
exploiting the fact that the length scale for the electric structure is 
very small compared to the length scale for magnetic structure. 

Finally, if the downstream state is Class D, there are two divergent 
states downstream, but only one free parameter upstream, so that 
longitudinal resistivity does not allow a shock to form in this case. In 
the next subsection we will discuss the meaning of the singular behavior in 
this regime. 








Calculating the Double Structured Shock 


As we have seen in the previous two subsections, when the plasma is 
making the transition from Class 8 to Class C, the electrostatic 
approximation breaks down and large electric fields and charge separations 
are set up by the flow. If the plasma were in Class C, and there were no 
dissipation, this charge separation results in an undamped ion acoustic 
wave. If it were in Class B with no dissipation the charge separation 
manifests itself as an ion acoustic soliton. With dissipation in the form 
of longitudinal resistivity, an ion acoustic shock can form so that the 
plasma can make a transition from Class B to Class C. Thus near this 
point, a narrow ion acoustic shock is embedded in the broader magnetic 
structure. This then corresponds to the isomagnetic potential jump 
discussed in Ref. 23. 

The choice of the electrostatic wave amplitude upstream is then 
equivalent to a choice of the position of the ion acoustic subshock within 
the broader magnetic structure. As we have seen, for a Class C plasma 
there is no magnetic perturbation which converges into the downstream side 

2 2 

of the ion acoustic subshock. This means for the case that k n >> u> V/vc , 

u pe 

where the magnetic structure and ion acoustic wave decouple, the ion 
acoustic wave must be at the position where the magnetic field has its 
downstream value. 

Hence the following scheme for computing the shock structure suggests 
itself. First integrate the equations for a quasi-neutral resistive shock, 
as written in Section III.A until one reaches the position where the 
magnetic field has the downstream value predicted by the Rankine-Hugoniot 
conditions. We find that at this point, the densities and temperatures are 
always below the downstream values, and also MV^ > 5/3 T g + y T i* At this 






point simply set the magnetic field equal to its constant value, set V ey * 
0 and solve for the ion acoustic shock in the unmagnetized plasma. The 
equations of the ion-acoustic shock are 


h (r "X + TT -iVi) ■ ”iV E + 14^ - v e) 2 - < 29 > 


The ion density and momentum equation are as in Eqs. (18) and (19). For 
the electrons 


Ve ‘ n oV (30 > 

fcV. --V E + “i-'’»(V 1 - T .h < 31 > 

h T n e V e T e * * ”e V e E + ITT “W V 1 ‘ < 32 > 

V 

and Poisson's Equation is still as given in Eq. (26). At the onset of the 
ion acoustic shock, the electric field can be initialized by taking its 
value from the quasi-neutral solution at the switch point. Of course, the 
exact value that E is initialized at is not very important because E 
reaches much higher values in the center of the transition region. 

Then Eqs. (18, 19, 26, 29-32) are integrated from the start of the ion 
acoustic shock on through to the downstream state. It is easy to show that 
this procedure preserves the overall Rankine-Hugoniot conditions; this is 
shown in Appendix A. 

To proceed, let us examine the singular nature of these equations by 

dT i 

writing out equations for V^. Eliminating from the energy equation 


(Eq. (29)) and inserting it in the momentum equation, Eq. (19), we find 








dx 


eV iE - mv jl (V i - V e )[ YVi - ( Y -l) 

;»i - ytj 


( V 1 - ».)! 


(33) 


Now the only singularity is at * yT^/M. It is especially easy to draw 
insight from this for the case Y ■ 3, one dimensional ions. As shown in 
Ref. 10, the fluid equations are those derived from a waterbag ion 
distribution for which the thermal width is (ST^/M)*^. Then the 
singularity at * 31^/M simply means that an ion at the outer edge of the 
waterbag distribution function has been reflected (its velocity has gone to 
zero in the shock frame). The obvious generalization then is to regard the 
2 

singularity at V - Y T^/M as the generalization of the condition for ion 
reflection to two or three dimensional ion distributions. Hence this 
singularity does have a simple interpretation in terms of reflected ions. 
In Ref. 10, it was shown that for ion acoustic shocks, one could add a 
reflected fluid upstream. The density of this reflected fluid then becomes 
an additional upstream parameter, which can be varied so as to eliminate 
another divergent downstream solution. This might also be possible to do 
for the magnetosonic shocks which we consider, but it will not be examined 
further here. 

To summarize, we have shown that longitudinal resistivity can give 
rise to double structured shocks in an intermediate range of Mach number. 
These are where the downstream state is Class C, as defined by Table 1. At 
the upper limit of this range, the shock structure once again becomes 
singular, but this singularity can now be interpreted in terms of reflected 
ions. 
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IV. Nonlocal Dissipation and Estimates for and a^. 

In the previous section, we discussed how longitudinal resistivity 
forms the transverse shock in the intermediate Mach number regime (Class 
C). The idea is that as the Mach number increases, large charge separation 
is generated, so an ion acoustic wave is set up. As this wave damps by 
longitudinal resistivity, the shock forms, as discussed in Section II. 
However, longitudinal resistivity is only a fluid model for the 
dissipation, the actual dissipation mechanism is electron or ion Landau 
damping, which is inherently non-fluid like. Nevertheless, it is an 
extremely powerful dissipative mechanism. If the ions are cold the damping 
distance due to electron Landau damping is extremely small, of order 
(M/m)*/ 2 kp 1 . This implies a of order If the ions are warm, the 
Landau damping rate is larger, implying still larger values of v^. 

As the longitudinal resistivity appears first in the momentum 
equation, we proceed by calculating the term in the electron momentum 
equation which arises from electron Landau damping. This means getting the 


correction to the pressure tensor. To do so we assume the plasma is 
infinite and homogeneous work in the (w, k) domain, and then write a 
linearized fluid equation. The only term which is not completely 
straightforward is the pressure tensor. As usual, we denote perturbed 
quantities with a superscript tilda. Then 


f 

e 


! 1 « 

eE 3v 
m -iui + ikv 


(34) 
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M -iui + ikv ’ 


where now a v denotes a particle velocity in a Vlasov model and, as usual, 
a capital V is the fluid velocity. The total perturbed electron density is 


■ / _ 
m J -iu) 


3v eEn 


[1 + /: 




(36) 


where we have assumed Maxwellian electrons, f g - (/2ir u @ ) 1 exp-(v 2 /2u^) and 
we work now in the plasma rest frame. Taking the resonant part of the 
second term in the brackets of Eq. (36), that is 


Im — + kV " Tri6(<j)~kv), 


(37) 


we find 


eE ri cir->l/2 

r— + nmn* 


(38) 


To continue, we would like to find the perturbed part of the total 
electron momentum flux 


f e - m / dvv 2 / dv v 3 k - • 


(39) 


The resonant part is small by a factor of (^—and we neglect it, so that 









To obtain an equation of state, we need the perturbed pressure in terms of 
the perturbed density. From Eqs. (38) and (40) we find 



(41) 


The force on the electrons is -ik?^. The first term in the square brackets 
3 >» 

simply gives - T g — n e> the pressure gradient force on an isothermal 
plasma. The second term, which is out of phase by n/2 is the additional 
dissipative contribution to the electron force. Using the fact that the 
linearized density equation gives the result n^ g “ u/k n e , we find the 
additional force density on the electrons is 

- ,2 

f - - ikP - - ikT n - n m u /*■ tut ^ • (42) 

e e e e e 2 |k| e v 

Modeling this additional retarding force as a longitudinal resistivity in 
the subshock, we find 


where we have assumed that within the subshock k ~ kp. 

Although the dissipation from Landau damping is quite simple in the k 
domain, it is very complicated and nonlocal in the x domain. In the k 

domain, the dissipative part of the perturbed pressure is the product of 
it k v 

-i nm u g | k| and V^. In the x domain, it is the convolution of the 
Fourier transforms of the two factors. Thus the dissipative part of the 
perturbed pressure in the x domain is 








V(x-) 


( 44 ) 


where P denotes the principal value.^ Thus the actual dissipation is 
nonlocal and quite complicated. 

Mow we consider the dissipative force on the ions. An exact analog of 
the calculation leading up to Eq. (42) gives the result that the 
dissipative force density on the ion fluid is given by 

< 45 > 

5f i 

For the case of positive phase velocity, < 0 for a thermal ion 

distribution, so that the force is still a retarding force. Assuming 
w/k is roughly constant as a function of k, the perturbed ion pressure in 
the x domain has the same form as given in Eq. (45). Depending on the ion 
temperature, the corresponding to the ion Landau damping can be large or 
small. If the electron temperature is comparable to the ion temperature, 
the damping length is comparable to k^, implying that 

Hence for either electron or ion dissipation, values for are quite 
large, easily large enough that the electrostatic potential is highly 
localized in the overall magnetic structure. Furthermore, the dissipation 
is not anomalous, in the sense that it does not arise from any instability, 
as does the dissipation which generates the transverse v. 

Since ion Landau damping heats the ions, and electron Landau damping 
heats the electrons, a reasonable guess for is 
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( 46 ) 


|ii 

_3v|v»V 


where ve have taken u/k - V, the flow speed. 

We now continue this section by discussing other aspects of the Landau 
damping model. As complicated as say Eq. (44) for P gd is, It is still a 
fairly coarse approximation to the actual nonlocal dissipation for several 
reasons. First of all, Eq. (44) is derived from linear theory, whereas the 
shock wave is inherently nonlinear. Secondly, since it is nonlocal, 

perturbed velocities everywhere contribute to the perturbed dissipative 
pressure at each point. However, only a portion of the plasma is 
electrostatic; a large part of it is dominated by the magnetic structure, 
and this was not treated correctly at all in the calculation of P gd . 
Third, the Fourier transform technique used to calculate P g works only 
if P g is either periodic or square integrable. However, a shock profile is 
neither so it is not clear just how valid the Fourier representation is. 
For instance, if V g (x) in Eq. (44) has a typical shock structure, then the 
integral in Eq. (44) diverges for large x. Thus the linear theory in the 
Fourier domain gives an infinite perturbed pressure at large x, that is, 
the region of x furthest from the shock. It is the divergence which is at 
the root of the claim by Pfirsch and Sudan^ that Landau damping cannot 
form a steady state shock. We feel that this divergence is clearly 
unphysical, so either nonlinear contributions play an important role in 
determining P e or else the Fourier transform scheme must be modified. 

Hence, since the actual dissipation mechanism is so complicated, a 
simple fluid approximation to it could be extremely useful. The only four 
possible conventional dissipation mechanisms are thermal conductivity, 
viscosity, transverse resistivity, and longitudinal resistivity. Thermal 








conduction can almost certainly be eliminated as a mechanism to form the 
shock because there does not seem to be any way that diffusing temperature 
can slow the forward motion and turn the upstream kinetic energy into 
downstream temperature. Viscosity can form the shock. However, there is 
no maximum Mach number for a viscous shock, no way that it is related to 
the breakdown of quasi-neutrality, and no way that it is related to ion 
reflection. Since all of these are known to be Important in experiments, 
viscosity does not appear to be the best fluid description of the 
dissipation either. Transverse resistivity is not only a good 
approximation, in all likelihood it is the dissipation mechanism, since 
transverse electron-ion streaming can give rise to instabilities which 
exchange momentum between electrons and ions. 7 However, as we have seen, 
transverse resistivity can only form the shock for Class B plasmas. 

This leaves longitudinal resistivity. While the actual dissipation 
mechanism is considerably more complex than this, it is a very good 
approximation to the physics. It is directly related to the breakdown of 
quasi-neutrality, it has an upper Mach number limit corresponding to ion 
reflection, it predicts a very narrow isomagnetic potential jump at the 
shock crest, and it has a lower Mach number limit corresponding to the 
maximum Mach number for transverse resistivity. All of these features are 
confirmed by the laboratory experiments we have cited. Furthermore, while 
longitudinal resistivity is simple, the two parameters which characterize 
it, and can be estimated from a coarse approximation to the actual 
dissipation mechanism. 

Let us conclude this section by discussing an important additional 
constraint on the physics. As we have seen, our basic model is that the 
flow steepens up until an ion acoustic wave is generated. The Landau 
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damping of this ion acoustic wave then provides the dissipation which forms 
the shock. However, for the theory of Landau damping to be valid ,^ the 
trapping time (for say the ions) 




1/2 


(47) 


where $ is the electrostatic potential, must be longer than the correlation 
time for a particle with velocity v. 


Z |E(k)| 2 «(kv-u>) 
k_ 

Z lE(k)I 2 
k 


(48) 


However, if the shock is simply a steady state potential profile, t = 00 

c 

for particles at the shock velocity. Hence the entire concept of Landau 
damping breaks down and the ion dynamics is dominated by trapping in a 
stationary potential. Thus in order for the Landau damping model to be 
valid, either the shock must have a temporally fluctuating structure, or 
else the total time of the experiment must be less than a trapping time. 
It is very likely that the shocks do have a fluctuating structure if only 
because in the broader part, dominated by transverse resistivity, the 
resistivity is probably caused by Debye length scale streaming 
instabilities. The spectrum of unstable waves will almost certainly 
overlap the electrostatic subshock and give rise to overall fluctuations in 


time 









V. Results 


In this section we present results of the calculations outlined in 
Section III. There are two classes of calculations. First, we consider 
parameters like those examined experimentally in Ref. 22. Here we model 
the longitudinal resistive heating as going entirely into the ions, as 
would be the case for strong ion Landau damping (that is not much less 
than T e ). Secondly, we consider a nitrogen plasma with parameters like 
that given in Ref. 21. More important, however, we consider no resistive 
ion heating. 

A. The Case of Strong Ion Heating 

In this case, we consider a shock in a hydrogen plasma with parameters 

as given in Ref. 22. That is, the upstream density is 4.6x10^ cm"-* and 

the upstream field is 730G. The ion y is taken as equal to 5/3, and T e ■ 

Tj - 1 ev. As is conventional in resistive shocks we assume a length of 

roughly 10 c/u) so that cv/w V. - 10 where we use upstream values for the 
pe pe A 

plasma frequency and Alfven speed. Although it is often assumed that for 
transverse resistive shocks there is no resistive ion heating, it is not 
easy to envision a transverse streaming instability which gives none. 
Also, a small amount of resistive ion heating would be extremely difficult 
to measure, and it might in fact be very important in determining the 
resistive shock width.® In addition, the ion temperature greatly affects 
the Ion Landau damping rate in the electrostatic portion of the shock. 
Accordingly we take a * 0.2, so there is 5 times as much resistive electron 
heating as resistive ion heating in the transverse part of the shock. 
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Since the transverse resistive shock always has a monatonic structure, 
the critical Mach number for transition to the Class C downstream state 
can be determined from the Rankine-Hugoniot conditions. For the 
parameters cited, this critical Mach number turns out to be 2.64. (Here 
the Mach number is defined as the shock velocity divided by the upstream 
magnetosonic speed). For Mach numbers below this, the shock is quasi¬ 
neutral and its structure is determined as specified in Section III.A. 

Above this Mach number, the quasi-neutral assumption breaks down. For 
larger Mach numbers, the shock structure is calculated as described in 
Secs. III.B and C. We assume, as seems consistent with experimental data, 
that in the electrostatic portion, ion Landau damping dominates. Thus we 
take * <*> and ■ y The upper critical Mach number, the 
transition from Class C to Class D cannot be determined from Rankine- 
Hugoniot conditions alone for two reasons. First of all the electrostatic 
part of the shock generally is not monatonic, but overshoots the final 
density, and then decays down to the downstream value. Secondly, even if 
the shock were monatonic, the Rankine-Hugoniot condition gives only T g + 
T^, whereas the transition from Class C to Class D depends only on T^. 
Therefore this upper critical Mach number depends on how the resistive 
heating is partitioned between electrons and ions. For the parameters we 
have chosen, we find that this upper critical Mach number turns out to be 
about 4.46. This agrees reasonably well with the upper limit Mach number 
of 4.25 +0.2 in Ref. 22. In Figs, la and lb are shown magnetic field and 
density profiles for a Class B shock with m ■ 2.5. Distances are given in 
centimeters. In Fig. 2a and 2b are shown magnetic field and density 
profiles in the transverse part of a supercritical shock having w * 4.4. 
The transverse part ends when the magnetic field achieves its downstream 
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value. At this point the electrostatic subshock is excited. The electron 
and ion density as a function of x in this region is shown in Fig. 2c. 
Again, distances are in centimeters. In all calculations, downstream 
values are carefully compared with Rankine-Hugoniot conditions as an 
independent check on the calculations. 

The actual ion heating downstream was measured in Ref. 22 by measuring 
n e (T g + T^) magnetically, and then measuring n £ and T g by laser 

scattering. The results, along with our calculated results, are shown in 
Fig. 3. The two basically agree, although the experiment shows somewhat 
more ion heating than what is calculated. 

Another important experimental result from Ref. 23 is that as the Mach 
2 2 

number increases, the value of 2e<|>/M(V o - Vj) decreases. Here <|> is the 
total potential jump and V Q and V d are upstream and downstream flow 
speeds. This goes hand in hand with the idea that as the Mach number 
increases, ions are slowed down not only by the electric field, but also by 
some additional dissipative effect. Our dissipative model, longitudinal 
resistivity also provides additional ion slowing down. In Fig. 4 are shown 
variations of this parameter with Mach number as measured by Eselevich et 
al., and as we calculated assuming upstream parameters in Ref. 22. Again 
the agreement is quite good. 

B. Shocks with No Ion Heating 

In this section we consider shocks with no ion heating, a * = 0. 

This corresponds to zero temperature ions, or else a waterbag ion 
distribution in one dimension. We consider the latter, so y^ ■ 3. Since 
the electrons and ions have different values of y, the Rankine-Hugoniot 
conditions cannot be calculated in the conventional way. However, by 
making use of the fact that u ■ 0, and the ions only heat adiabatically, 


35 







Che conservation relations can be calculated. This is done in Appendix 
C. For upstream fluid parameters, we use those in Ref. 21, a singly 
ionized nitrogen plasma, n * lO 1 ^, B = 2KG, T e = 2 ev, T i = 0.5 ev, 
cv/uip e V A * 10, and =■ ^pe’ corres P onc **- n 8 to electron-Landau damping. 
Then, Rankine-Hugoniot conditions show that the Mach number for transition 
from Class B to Class C is M = 2.68. There is no second critical Mach 
number for transition from Class C to Class D predited by Rankine-Hugoniot 
simply because there is not enough ion heating. However, the shock is not 
monotonic, but oscillatory, so that in an overshoot it is possible to reach 
the second critical point * 3T^/M. We find this singular point is 
reached at m * 3.73. In Fig. 5a and 5b are shown magnetic field and 
density profiles versus position in the quasi-neutral part for a shock 
with m m 3.5. In Fig. 5c is the ion density versus x in the electrostatic 






VI 


Conclusions 


Our first conclusion is that the experimental evidence strongly 
indicates the existence of a significant intermediate range of Mach numbers 
between shocks formed by transverse resistivity and shocks formed by 
reflected ions. In this range, according to Ref. 23, the breakdown of 
quasi-neutrality plays an important role and an electrostatic subshock 
forms on the crest of the magnetic structure. Second, we have shown that 
the singular behavior of the quasi-neutral fluid equations is not related 
to ion reflection but rather to the breakdown of quasi-neutrality. In the 
intermediate regime then, an ion acoustic wave train is set up. We 
speculate that the Landau damping of this acoustic wave is the dissipation 
mechanism which forms the shock. Third, we have investigated longitudinal 
resistivity as a simple fluid model for the dissipation mechanism for shock 
formation. We conclude that it is a reasonable model. It is important 
only between the upper Mach number for transverse resistivity, and the 
lower Mach number for a fluid picture of ion reflection. It also is 
directly related to the breakdown of the quasi-neutrality assumption, and 
it predicts an electrostatic subshock at the crest of the magnetic 
structure. 
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Figure 2 


a) Magnetic field and, b) density profiles in the trans\ 
part of the supercritical shock having m ■ 4.4. c) 
(solid) and electron (dotted) density profile in the elec 


static part of the shock. 













Figure 3 Calculated and measured downstream ion temperatures 


intermediate range shocks as a function of m . 



























Appendix A 


In this appendix, we show that the scheme for determining the shock 
structure, as specified in Section III.C satisfy the overall conservation 
relations. Three of the conservation relations, that for mass, momentum 
and energy are satisfied at each point in the shock. The fourth relation 
uB - U q B 0 is only satisfied downstream. 

Consider now the conservation equations for momentum and energy 
flux. Each flux is the sum of a fluid and electromagnetic component. In 
the quasi-neutral region, the electrostatic part of the momentum flux is 
neglected, but the magnetic part of the momentum flux, B^/8 it and the 
electromagnetic part of the energy flux ^ BE q are each included in the 
formulation. In the quasi-neutral region then, n, V, and T are determined, 
in part, so that particle, momentum and total energy flux are conserved. 
At the switch point, B takes on its downstream value and remains 
constant. Therefore in the electrostatic region, the magnetic part of the 
momentum flux and electromagnetic part of the energy flux are constant and 
have the proper downstream value. 

In the electrostatic region, the formulation leads to a conservation 

law for total energy flux and a conservation law for total fluid momentum 

E2 

flux plus electrostatic momentum flux - -g^- . Downstream, of course, E • 0 
and the fluid momentum and energy fluxes are as given by the Rankine- 
Hugoniot condition. 

Now consider the fluid energy flux as one moves back upstream toward 
the switch point in the electrostatic region. Since it is constant, it 
matches smoothly onto that in the quasi-neutral regime. Next consider the 
momentum flux. Since the formulation for the electrostatic region contains 
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a momentum flux - -g^- , and the quasi-neutral formulation does not, the 

momentum fluxes do not quite natch at the switch point if they have the 

E 2 

same values downstream. However, if -g— is a very small part of the 

momentum flux at the switch point, the error in momentum flux is very 

small. Using the expression for E at the switch point, Eq. (2), and 

considering say the ratio of electrostatic to thermal momentum flux, we 

find /nT >. (u? v/vc^k « 1. In all of our calculations of shock 

8it e v pe D- 1 

structure, we have checked that the electrostatic component of the momentum 
flux is negligible at the switch point, so that total momentum is conserved 
throughout the shock. Of course, while the electrostatic momentum flux is 
very small at the beginning and end of the electrostatic shock, it can be 
very large at points in between. 
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In this appendix, we calculate Rankine-Hugoniot conditions for the 


case of ■ y , « 3, but where the ions are only heated adiabatically 

so that ■ T^ 0 Then denoting upsteam quantities with a subscript 

1 o 

zero and downstream quantities with no subscript, the conservation 
equations for momentum and energy flux reduce to two equations for n and 

V 



io 3 O/n ^2 

~T~ n + 17 t") 


I (—) 2 
io '•n 1 


' 2 A io 


Bn 

o 

4irn^ 


(B2) 


In writing Eqs. (Bl) and (B2) we have used the fact that n Q V 0 ■ nV and V Q B 0 
• VB. Now each equation is linear in T @ , so solving for T g from one and 
inserting in the other, we find 


oi , MV* i 

•"ft 1 - b + W 1 -») + -r 2 - ( l - “) + ( l - ») T e[ > ‘ 

MV 2 

| “hr t 1 - -r) + t V 1 - » 2 ) + "'lot 1 - “)1 < B3 > 

u 


where u ■ — and V, * B /4irn M. Factoring out the obvious root u * 1, we 
n Ao o o ° 

o 

arrive at a cubic equation for u 



(B4) 












This equation also has a root u * 1 if MV * MV, + T T + 3T. as can be 

o Ao J eo i 

easily verified. The procedure then is to solve Eq. (B4) numerically 
starting from this root and following it as V Q increases. Assuming a 
singly ionized nitrogen plasma with n Q - 10 15 , T go = 2ev, T io = 0.5ev, B * 
2KG, Fig. B1 shows a plot of downstream values of density and electron 
temperature versus magnetosonic Mach number m - v /fV? + 4 T /M + 









(solid) for a shock with y g * 5/3, y. 
heating. 


3 and adiabatic ion 
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220 7TH STREET, NE 
CHARLOTTESVILLE, VA 22901 
01CY ATTN DRXST-SD 

COMMANDER 

U.S. ARMY MATERIAL DEV & READINESS CMD 
5001 EISENHOWER AVENUE 
ALEXANDRIA, VA 22333 

01CY ATTN DRCLDC J.A. BENDER 

COMMANDER 

U.S. ARMY NUCLEAR AND CHEMICAL AGENCY 
7500 BACKLICK ROAD 
BLDG 2073 

SPRINGFIELD, VA 22150 
01CY ATTN LIBRARY 


DIRECTOR 

U.S. ARMY BALLISTIC RESEARCH LABORATORY 
ABERDEEN PROVING GROUND, MD 21005 

01CY ATTN TECH LIBRARY EDWARD EAICY 

COMMANDER 

U.S. ARMY SATCOM AGENCY 
FT. MONMOUTH, NJ 07703 

01CY ATTN DOCUMENT CONTROL 

COMMANDER 

U.S. ARMY MISSILE INTELLIGENCE AGENCY 
REDSTONE ARSENAL, AL 35809 
01CY ATTN JIM GAMBLE 

DIRECTOR 

U.S. ARMY TRADOC SYSTEMS ANALYSIS ACTIVITY 
WHITE SANDS MISSILE RANGE, NM 88002 
01CY ATTN ATAA-SA 
01CY ATTN TCC/F. PAYAN JR. 

01CY ATTN ATTA-TAC LTC J. HESSE 

COMMANDER 

NAVAL ELECTRONIC SYSTEMS COMMAND 
WASHINGTON, D.C. 20360 

01CY ATTN NAVALEX 034 T. HUGHES 
01CY ATTN PME 117 
01CY ATTN PME 117-T 
01CY ATTN CODE 5011 

COMMANDING OFFICER 
NAVAL INTELLIGENCE SUPPORT CTR 
4301 SUITLAND ROAD, BLDG. 5 
WASHINGTON, D.C. 20390 

01CY ATTN MR. DUBBIN STIC 12 

01CY ATTN NISC-50 

01CY ATTN CODE 5404 J. GALET 

COMMANDER 

NAVAL OCCEAN SYSTEMS CENTER 
SAN DIEGO, CA 92152 
01CY ATTN J. FERGUSON 







NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 

01CY ATTN CODE 4700 S. L. Ossakou 

26 CYS IF UNCLASS. I CY IF CLASS) 
01CY ATTN CODE 4701 I Vitkovitsky 

01CY ATTN CODE 4780 J. Huba (100 

CYS IF UNCLASS, 1 CY IF CLASS) 
01CY ATTN CODE 7500 

01CY ATTN CODE 7550 

01CY ATTN CODE 7580 

01CY ATTN CODE 7551 

01CY ATTN CODE 7555 

01CY ATTN CODE 4730 E. MCLEAN 

01CY ATTN CODE 4108 

01CY ATTN CODE 4730 B. RIPIN 

20CY ATTN CODE 2628 

COMMANDER 

NAVAL SEA SYSTEMS COMMAND 
WASHINGTON, D.C. 20362 
01CY ATTN CAPT R. PITKIN 


COMMANDER 

AEROSPACE DEFENSE COMMAND/DC 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 
01CY ATTN DC MR. LONG 

COMMANDER 

AEROSPACE DEFENSE COMMAND/XPD 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 
OICY ATTN XPDQQ 
01CY ATTN XP 


AIR FORCE GEOPHYSICS LABORATORY 
RANSCOM AFB, MA 01731 

OICY ATTN OPR HAROLD GARDNER 
OICY ATTN LKB KENNETH S.W. CHAMPION 
OICY ATTN OPR ALVA T. STAIR 
OICY ATTN PHD JURGEN BUCHAU 
OICY ATTN PHD JOHN P. MULLEN 


COMMANDER 

NAVAL SPACE SURVEILLANCE SYSTEM 
DAHLGREN, VA 22448 

OICY ATTN CAPT J.H. BURTON 

OFFICER-IN-CHARGE 
NAVAL SURFACE WEAPONS CENTER 
WHITE OAK, SILVER SPRING, MD 20910 
OICY .ATTN CODE F31 

DIRECTOR 

STRATEGIC SYSTEMS PROJECT OFFICE 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20376 
OICY ATTN NSP-2141 
OICY ATTN NSSP-2722 FRED WIMBERLY 

COMMANDER 

NAVAL SURFACE WEAPONS CENTER 
DAHLGREN LABORATORY 
DAHLGREN, VA 22448 

OICY ATTN CODE DF-14 R. BUTLER 

OFFICER OF NAVAL RESEARCH 
ARLINGTON, VA 22217 
OICY ATTN CODE 465 
OICY ATTN CODE 461 
OICY ATTN CODE 402 
OICY ATTN CODE 420 
OICY ATTN CODE 421 


AF WEAPONS LABORATORY 
KIRTLAND AFT, NM 87117 
OICY ATTN SUL 

OICY ATTN CA ARTHUR H. GUENTHER 
OICY ATTN NTYCE 1LT. G. KRAJEI 

AFTAC 

PATRICK AFB, FL 32925 
OICY ATTN TF/MAJ WILEY 
OICY ATTN TN 

AIR FORCE AVIONICS LABORATORY 
WRIGHT-PATTERSON AFB, OH 45433 
OICY ATTN AAD WADE HUNT 
OICY ATTN AAD ALLEN JOHNSON 

DEPUTY CHIEF OF STAFF 
RESEARCH, DEVELOPMENT, & ACQ 
DEPARTMENT OF THE AIR FORCE 
WASHINGTON, D.C. 20330 
OICY ATTN AFRDQ 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 
OICY ATTN J. DEAS 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/YSEA 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01732 




HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/DC 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 

OICY ATTN DCKC MAJ J.C. CLARK 

COMMANDER 

FOREIGN TECHNOLOGY DIVISION, AFSC 
WRIGHT-PATTERSON AFB, OH 45433 
01CY ATTN NICD LIBRARY 

01CY ATTN ETDP B. BALLARD 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
GRIFFISS AFB, NY 13441 

01CY ATTN DOC LIBRARY/TSLD 
01CY ATTN OCSE V. COYNE 

SAMSO/SZ 

POST OFFICE BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
(SPACE DEFENSE SYSTEMS) 

01CY ATTN SZJ 

STRATEGIC AIR COMMAND/XPFS 
OFFUTT AFB, NB 68113 

01CY ATTN ADWATE MAJ BRUCE BAUER 
01CY ATTN NRT 

01CY ATTN DOK CHIEF SCIENTIST 

SAMSO/SK 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 

01CY ATTN SKA (SPACE COMM SYSTEMS) 
M. CLAVIN 

SAMSO/MN 

NORTON AFB, CA 92409 
(MINUTEMAN) 

01CY ATTN MNNL 


COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
HANSCOM AFB, MA 01731 

01CY ATTN EEP A. LORENTZEN 


DEPARTMENT OF ENERGY 
LIBRARY ROOM C-042 
WASHINGTON, D.C. 20545 

01CY ATTN DOC CON FOR A. LABOWITZ 


DEPARTMENT OF ENERGY 
ALBUQUERQUE OPERATIONS OFFICE 
P.O. BOX 5400 
ALBUQUERQUE, NM 87115 

01CY ATTN DOC CON FOR D. SHERWOOD 

ECS.G, INC. 

LOS ALAMOS DIVISION 

P.O. BOX 809 

LOS ALAMOS, NM 85544 

01CY ATTN DOC CON FOR J. BREEDLOVE 


UNIVERSITY OF CALIFORNIA 
LAWRENCE LIVERMORE LABORATORY 
P.O. ..OX 808 
LIVERMORE, CA 94550 

01CY ATTN DOC CON FOR TECH INFO DEPT 
01CY ATTN DOC CON FOR L-389 R. OTT 
01CY ATTN DOC CON FOR L-31 R. HAGER 
01CY ATTN DOC CON FOR L-46 F. SEWARD 

LOS ALAMOS NATIONAL LABORATORY 

P.O. BOX 1663 

LOS ALAMOS, NM 87545 

01CY ATTN DOC CON FOR J. WOLCOTT 
01CY ATTN DOC CON FOR R.F. TASCHEK 
OICY ATTN DOC CON FOR E. JONES 
01CY ATTN DOC CON FOR J. MALIK 
OICY ATTN DOC CON FOR R. JEFFRIES 
OICY ATTN DOC CON FOR J. ZINN 
OICY ATTN DOC CON FOR P. KEATON 
OICY ATTN DOC CON FOR D. WESTERVELT 
OICY ATTN D. SAPPENFIELD 

SANDIA LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 

OICY ATTN DOC CON FOR W. BROWN 

OICY ATTN DOC CON FOR A. THORNBROUCH 

OICY ATTN DOC CON FOR T. WRIGHT 

OICY ATTN DOC CON FOR D. DAHLGREN 

OICY ATTN DOC CON FOR 3141 

OICY ATTN DOC CON FOR SPACE PROJECT DI? 

SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE,- CA 94550 

OICY ATTN DOC CON FOR B. MURPHEY 
OICY ATTN DOC CON FOR T. COOK 

OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON DR. YO SONC 
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OTHER GOVERNMENT 


BOEING COMPANY, THE 
P.O. BOX 3707 

DEPARTMENT OF COMMERCE SEATTLE, UA 98124 

NATIONAL BUREAU OF STANDARDS 01CY ATTN C. KEISTER 

WASHINGTON, D.C. 20234 OICY ATTN D. MURRAY 

01CY (ALL CORRES: ATTN SEC OFFICER FOR) OICY ATTN G. HALL 

OICY ATTN J. KENNEY 


INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & INFO ADMIN 
BOULDER, CO 80303 

OICY ATTN A. JEAN (UNCLASS ONLY) 

OICY ATTN W. UTLAUT 
OICY ATTN D. CROMBIE 
OICY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
OICY ATTN R. GRUBB 
OICY ATTN AERONOMY LAB G. REID 


DEPARTMENT OF DEFENSE CONTRACTORS 


AEROSPACE CORPORATION 
P.O. BOX 92957 
LOS ANGELES, CA 90009 
OICY ATTN I. GARFUNKEL 
OICY -ATTN T. SALMI 
OICY ATTN V. JOSEPHSON 
OICY ATTN S. BOWER 
OICY ATTN D. OLSEN 

ANALYTICAL SYSTEMS ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 

OICY ATTN RADIO SCIENCES 

AUSTIN RESEARCH ASSOC., INC. 

1901 RUTLAND DRIVE 
AUSTIN, TX 78758 
OICY ATTN L. SLOAN 
OICY ATTN R. THOMPSON 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 
OICY ATTN J. WORKMAN 
OICY ATTN C. PRETTIE 
OICY ATTN S. BRECHT 


CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
OICY ATTN D.B. COX 
OICY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LISTHICUM ROAD 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 


CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 

OICY ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
OICY ATTN H. LOCSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

EOS TECHNOLOGIES, INC. 

606 Wllshire Blvd. 

Santa Monica, Calif 90401 
OICY ATTN C.B. GABBARD 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M.H. BORTNER SPACE SCI LAB 

GENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSE, NY 13201 
OICY ATTN F. REIBERT 








GENERAL ELECTRIC TECH SERVICES CO., INC. 
HUES 

COURT STREET 
SYRACUSE, NY 13201 
01CY ATTN G. HILLMAN 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER) 
01CY ATTN T.N. DAVIS (UNCLASS ONLY) 
01CY ATTN TECHNICAL LIBRARY 
01CY ATTN NEAL BROWN (UNCLASS ONLY) 

GTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

01CY ATTN DICK STEINHOF 

HSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

01CY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

01CY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
1801 NO. BEAUREGARD STREET 
ALEXANDRIA, VA 22311 
Q1CY ATTN J.M. AEIN 
01CY ATTN ERNEST BAUER 
01CY ATTN HANS WOLFARD 
01CY ATTN JOEL BENGSTON 


INTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 

01CY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 
P.O. BOX 85154 
SAN DIEGO, CA 92138 

01CY ATTN J.L. SPERLING 


JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 

01CY ATTN DOCUMENT LIBRARIAN 
01CY ATTN THOMAS POTEMRA 
OICY ATTN JOHN DASSOULAS 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN T. MEAGHER 

KAMAN TEMPO-CENTER FOR ADVANCED STUDIES 
816 STATE STREET (P.O DRAWER QQ) 

SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 

LINXABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 

LOCKHEED MISSILES & SPACE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES 6 SPACE CO., INC. 

3251 HANOVER STREET 
PALO ALTO, CA 94304 
OICY ATTN MARTIN WALT DEPT 52-12 
OICY ATTN W.L. IMHOF DEPT 52-12 
OICY ATTN RICHARD G. JOHNSON DEPT 52-12 
OICY ATTN J.B. CLADIS DEPT 52-12 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, FL 32805 
OICY ATTN R. HEFFNER 


M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN L. LOUGHLIN 
OICY ATTN D. CLARK 








MCDONNEL DOUGLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
01C? ATTN N. HARRIS 
Old ATTN J. MOULE 
Old ATTN GEORGE MROZ 
Old ATTN W. OLSON 
Old ATTN R.W. HALPRIN 
Old ATTN TECHNICAL LIBRARY SERVICES 

MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
Old ATTN P. FISCHER 
Old ATTN W.F. CREVIER 
Old ATTN STEVEN L. GUTSCHE 
Old ATTN R. BOGUSCH 
Old ATTN R. HENDRICK 
Old ATTN RALPH KILB 
Old ATTN DAVE SOULE 
Old ATTN F. FAJEN 
Old ATTN M. SCHEIBE 
Old ATTN CONRAD L. LONGHIRE 
Old ATTN B. WHITE 

MISSION RESEARCH CORP. 

1720 RANDOLPH ROAD, S.E. 

ALBUQUERQUE, NEW MEXICO 87106 
01CY R. STELLINCWERF 
Old M. ALME 
Old L. WRIGHT 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

Old ATTN JOHN MORGANSTERN 
Old ATTN G. HARDING 
Old ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
01CY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH COR? 

12340 SANTA MONICA BLVD. 

LOS ANGELES, CA 90025 

01CY ATTN E.C. FIELD, JR. 


PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

01CY ATTN IONOSPHERIC RESEARCH LAB 

PHOTOMETRICS, INC. 

4 ARROW DRIVE 
WOBURN, MA 01801 

01CY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

Old ATTN E.J. FREMOUW 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

R 6 D ASSOCIATES 
P.O. BOX 9695 
MARINA DEL REY, CA 90291 
01CY ATTN FORREST GILMORE 
OICY ATTN WILLIAM B. WRICHT, JR. 
Old ATTN ROBERT F. LELEVIER 
OICY ATTN WILLIAM J. KARZAS 
OICY ATTN H. ORY 
OICY ATTN C. MACDONALD 
Old ATTN R. TURCO 
OICY ATTN L. DeRAND 
OICY ATTN W. TSAI 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
OICY ATTN CULLEN CRAIN 
OICY ATTN ED BEDROZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

OICY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
330 WEST 42nd STREET 
NEW YORK, NY 10036 

OICY ATTN VINCE TRAPANI 
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SCIENCE APPLICATIONS, INC. 

1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 

01CY ATTN LEWIS M. LINSON 
01CY ATTN DANIEL A. HAMLIN 
01CY ATTN E. PRIEMAN 
01CY ATTN E.A. STRAKER 
01CT ATTN CURTIS A. SMITH 
01CY ATTN JACK MCDOUCALL 

SCIENCE APPLICATIONS, INC 
1710 COODRIDGE DR. 

MCLEAN, VA 22102 
ATTN: J. COCKAYNE 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 

01CY ATTN DONALD NEILSON 
01CY ATTN ALAN BURNS 
01CY ATTN G. SMITH 
01CY ATTN R. TSUNODA 
01CY ATTN DAVID A. JOWSON 
01CY ATTN WALTER G. CHESNUT 
01CY ATTN CHARLES L. RINO 
01CY ATTN WALTER JAYE 
01CY ATTN J. VICKREY 
01CY ATTN RAY L. LEADABRAND 
01CY ATTN G. CARPENTER 
01CY ATTN G. PRICE 
01CY ATTN R. LIVINGSTON 
01CY ATTN V. GONZALES 
01CY ATTN D. MCDANIEL 

TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 

01CY ATTN W.P. BOQUIST 

TOYON RESEARCH CO. 

P.O. Box 6890 
SANTA BARBARA, CA 93111 
01CY ATTN JOHN ISE, JR. 

01CY ATTN JOEL GARBARINO 


TRW DEFENSE & SPACE SYS GROUP 
ONE SPACE PARK 
REDONDO BEACH, CA 90278 
01CY ATTN R. K. PLEBUCH 
01CY ATTN S. ALTSCHULER 
01CY ATTN D. DEE 
01CY ATTN D/ STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
01CY ATTN W. REIDY 
01CY ATTN J. CARPENTER 
01CY ATTN C. HUMPHREY 
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